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Abstract. Relativistic shocks are usually thought to occur in violent astrophysical explosions. 
These collisionless shocks are mediated by a plasma kinetic streaming instability, often loosely 
referred to as the Weibel instability, which generates strong magnetic fields "from scratch" very 
efficiently. In this review paper we discuss the shock micro-physics and present a recent model of 
"pre-conditioning" of an initially unmagnetized upstream region via the cosmic -ray-driven Weibel- 
type instability. 
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INTRODUCTION 

It has been shown in recent years that collisionless relativistic shocks are mediated by 
the Weibel instability — a current filamentation instability that produces strong, sub- 
equipartition magnetic fields at the shock front [1]. This is an attractive model for 
gamma-ray bursts (GRBs), because it puts a synchrotron shock model on a firm physical 
ground. Here we describe the physics of the instability, its linear regime and saturation. 
We discuss the properties of the nonlinear Weibel turbulence — the state of the ongoing 
self- similar process of current mergers. Finally, we present a toy model of the foreshock 
region — an extended region in the upstream of the shock which may likely be populated 
with strong and relatively large-scale (meso- scale) magnetic fields produced by shock- 
accelerated particles. The foreshock can be the main player in determining the shock 
radiative processes and the efficiency of cosmic ray acceleration. 

LINEAR REGIME OF FIELD GROWTH AND ITS SATURATION 

The instability under consideration was first predicted by Weibel (2D for a non-relativistic 
plasma with an anisotropic distribution function. The simple physical interpretation con- 
siders the PDF anisotropy more generally as a two-stream configuration of cold plasma. 
Below we give a brief, qualitative description of this streaming magnetic instability. 

Let us consider, for simplicity, the dynamics of one species only (e.g., protons), 
whereas the other (electrons) is assumed to provide global charge neutrality^ The 
streaming particles are assumed to move along the z-axis with the velocities v = +zv z 
and v = — zv z , thus forming equal particle fluxes in opposite directions (so that the net 



In reality, the role of protons is more complicated, e.g., they play a crucial role in the electron heating, 
which we do not consider here. 



current is zero). Such a particle distribution occurs naturally near the front of a shock 
(moving along z-direction), where the "incoming" (in the shock frame) ambient gas par- 
ticles meet the "outgoing" particles reflected from the shock potential (loosely speaking, 
the low energy cosmic rays). Thus, the particle velocities v z are of order the shock veloc- 
ity, v z ~ v s h- The counter- streaming particles may also have some thermal spread. Since 
for high-Mach number shocks, v t h <C v s h in the upstream region, we may neglect the par- 
allel velocity spread in our consideration. The thermal spread in the transverse direction 
cannot be neglected, however. We parameterize the PDF anisotropy as follows: 
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where £m oc v 2 ~ v 2 h is the energy of particle along z-direction, oc (y 2 + v 2 ) oc v thermal 
c 2 . is the thermal energy in the transverse direction, £ tot = ey + is the total energy, c s 
is the sound speed upstream and the Mach number of the shock is M = v s h/c y . Clearly, 
for strong shocks M» 1, the anisotropy parameter is close to unity, A ~ 1. Next, 
according to the linear stability analysis technique, we add an infinitesimal magnetic 
field fluctuation, B = xB x cos(ky). The Lorentz force, e{\ x B)/c, acts on the charged 
particles and deflects their trajectories, as is shown in Figure la. As a result, the protons 
moving upward and those moving downward will concentrate in spatially separated 
current filaments. The magnetic field of these filaments appears to increase the initial 
magnetic field fluctuation. The growth rate and the wavenumber of the fastest growing 
mode (which, in fact, sets the spatial correlation scale of the produced field) are 

y B = A (o P)S {v z /c) , k B =A (Dpj/c, (2) 

where co p , = (4xe 2 n s /m s ) ^ is the plasma frequency defined for species , (electrons, 
protons, etc.), n p and m p are the number density and the mass of the protons, respec- 
tively. The above scalings are for non-relativistic plasmas. In the relativistic case, one 
shall simply substitute m s — > y s m s , there y s is the Lorentz factor of the streaming species. 
Note that the instability is driven by the PDF anisotropy and should quench for the 
isotropic case, A = 0. 

The Lorentz force deflection of particle orbits increases as the magnetic field pertur- 
bation grows in amplitude. The amplified magnetic field is random in the plane per- 
pendicular to the particle motion, since it is generated from a random seed field. Thus, 
the Lorentz deflections result in a pitch angle scattering, which makes the bulk of the 
PDF isotropic. If one starts from a strong anisotropy, so that the thermal spread is much 
smaller than the particle bulk velocity, most of the particles will eventually isotropize 
and the thermal energy associated with their random motions will be equal to their ini- 
tial directed kinetic energy. This final state will bring the instability to saturation. 

The saturation level of the magnetic field may readily be estimated as follows. First 
of all, note that the instability is due to the free streaming of particles. As the magnitude 
of the magnetic field grows, the transverse deflection of particles gets stronger, and their 
free streaming across the field lines is suppressed, see Figure lb. Once the bounce rate 
of the streaming particles in the self-generated current filaments becomes greater than 
the instability growth rate, the instability ceases. In the relativistic case we are concerned 
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FIGURE 1. Illustration of various stages of the Weibel instability. Color coding of particles: blue — the 
incoming particles from the IGM, red — the particles scattered from the shock, (a) Linear regime: current 
filamentation; (b) saturation; (c) nonlinear regime: filament coalescence. 



with here, this condition is identical to that the typical curvature scale for the deflections 
is the Larmor radius, 

Pl = V^ B /G)c,s: (3) 

where v±b is the particle velocity transverse to the direction of the local magnetic field 
and 

eB 

(0c, s = (4) 

y s m s c 

is the cyclotron (Larmor) frequency of species s and we assumed v±b ~ c. On scales 
larger than pi, particles can only move along field lines. Hence, when the growing 
magnetic fields become such that kgpt ~ 1, the particles are magnetically trapped and 
can no longer amplify the field. Assuming an isotropic particle distribution at saturation 
(v±b ~ v s h), this condition can be re- written as 

£B = T" ^ ? ~ A ■ ^ 

m * C "Whack 

For strong shocks (M>1,A~1), this corresponds to the magnetic energy density close 
to equipartition with the thermal energy of particles downstream the shock. 

Numerical PIC simulation indicate a somewhat smaller value of typically not 
exceeding 0.1 or so even for strong shocks, that is the magnetic field energy density 
at a relativistic shock is of the order of 10% of the shock kinetic energy [@,M[5|]. For 
more recent relevant work, see [@,0]. 



NONLINEAR WEIBEL TURBULENCE AND FILAMENT 

COALESCENCE 

Numerous 3D PIC numerical simulations demonstrate that the generated magnetic fields 
are associated with a quasi-two-dimensional distribution of current filaments. Hence 
we suggest the following toy model. We consider straight one-dimensional current 
filaments oriented in the vertical, z-direction. Initially, all filaments are identical: the 
initial diameter of them is Do, their initial mass per unit length is /io ~ mn{%D^/A), 
where m is the mass of plasma particles (e.g., electrons) and n is their number density. 
Each filament carries current Iq in either positive or negative f-direction. The net current 
in the system is set to zero, i.e., there are equal numbers of positive and negative 
current filaments. We also assume that the initial separation (the distance between the 
centers) of the filaments do is comparable to their size, do — 2Dq. Finally, no external 
homogeneous magnetic field is present in the system. Here we assume that filaments 
have a simple structure: they have no return current on the outside. Generally, this is not 
true, as the filaments represent some sort of the Harris equilibrium. The return current 
can be incorporated in the analysis and is expected to decrease the overall coalescence 
rate, because it shall partially screen the main current and, thus, reduce the interaction 
strength between the filaments. 

Initially, the filaments are at rest and their positions in space are random. This con- 
figuration is unstable because opposite currents repel each other, whereas like currents 



are attracted to each other and tend to coalesce and form larger current filaments. The 
characteristic scale of the magnetic field will accordingly increase with time. We study 
this process quantitatively using the toy model of two interacting filaments. 

The magnetic field produced by a straight filament is Bo(r) = 2/o/(cr), where r 
is the cylindrical radius. The force per unit length acting on the second filament is 
dF/dl = —BqIq/c. Since dF jdl = /Ax, where x is the position in the center of mass 
frame and "overdot" denotes time derivative, we write the equation of motion as follows: 
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where we used that r = 2x and the reduced mass jl r = Hq/2. We define the coalescence 
time as the time required for the filaments, which are initially at rest, to cross the distance 
between them and "touch" each other, which happens when the distance between their 
centers becomes equal to Do, i.e., when x = Dq/2. The coalescence time, as it is defined 
above, is independent of the details of the merging process itself, which involves rather 
complicated dynamics associated with the redistribution of currents. Quite obviously, 
the interaction between the filaments is the weakest at large distances x ~ xq ~ Jo/2. 
Hence, the coalescence rate is limited by the filament motions at the largest scales. The 
coalescence time can be readily estimated from Eq. ©, assuming that x ~ xq ~ do/2 
and x ~ (do/2)r 2 , as follows: 

To,NR~{D 2 c 2 iio/(2lZ)) l/2 . (7) 

The above estimate is valid as long as the motion is non-relativistic. The maximum 
velocity of a filament is at the time of coalescence, x = Dq/2: 

v ~D /2t ~I /(c^/2^). (8) 

It must always be much smaller than the speed of light. 

If the motion of a filament during the merger becomes sub-relativistic, i.e., v merger 
becomes a significant fraction of c the separation cannot decrease faster than as t(x) ~ 
x/c. Therefore, the coalescence time will be 

T 0r R ^Dq/ V me rger- (9) 

The filament coalescence is a hierarchical process. Indeed, suppose that initially the 
system contains Nq current filaments, with an average separation do ~ 2Dq. Each of the 
filaments carries current Iq, its diameter is Do and its mass per unit length is /Jo- For 
simplicity, we assume that filaments coalesce pairwise. 

Having the original "zeroth generation" of filaments merged (the process takes about 
To,NR or ^o,R t0 complete), the system will now contain Nq/2 of "first generation" 
filaments. Each of these filaments carries current I{ = 2Iq, has mass per unit length 
Hi = 2/j.q, and the separation between them is d\ = V2do (because the two dimensional 
number density of filaments decreased by 2). Since ,u « D 2 , the filament size also 
increases as D\ = V2Dq. Remarkably, this new configuration is identical to the initial 



FIGURE 2. The 2D structure of the magnetic fields in e~e + plasmas at various times (from H]). 
The change of the field correlation length with time is clearly seen. Similar The growth of this length 
is substantially slower and the magnetic field filling factor is respectively larger in the electron-proton 
plasma. 



one, but with the re- scaled parameters. Hence, the coalescence process is self- similar. 
The produced first generation filaments will be interacting with each other and merge 
again to yield the second generation. The coalescence process will then continue in a 
self-similar way. Note that the coalescence times at each stage are not necessarily the 
same. Taking into account that at the k-ih. merger level, i.e., after k pairwise mergers: 
4 = 2%, ju* = 2*/io, D k = 2 k / 2 D , d k ~ D k /2, we obtain 

*k,NR = ?0,NR, %R = 2 k l 2 T^ R . (10) 

Since the coalescence time is independent of k while the filaments are non-relativistic, 
whereas the distance between them increases, the typical velocities of the merging 
filaments grow with time and, will approach some terminal velocity v merger ~ c, as 
v k = 2 k ' 2 v Q . 

Finally, it is instructive to present the evolution of the parameters as a function of 
physical time, t, rather than the merger level, k. Apparently, it takes t = E^ =0 % to 
complete k mergers, where % k is given by Eq. (flOl) . Thus, for the non-relativistic and 
relativistic cases respectively, we have: k = t/tQ^R and k ~ 21og 2 [(t/^ofi)]- Thus, the 
magnetic field correlation length increases as a function of time as 

X B (t) = Do^A 2 ^, X B ( t ) ~ v merger t, (11) 

in the non-relativistic and sub-relativistic regimes, respectively. Note that the last ex- 
pression is an approximation at large times t ^> Tq^r, i.e., at large k ^> 1. 

Simulations of the Weibel filament dynamics in 2D pair-dominated and electron- 
ion plasmas confirm the above arguments. A few snapshots are shown in Figure [2l It 
has been observed that the filament coalescence does occur and the merger velocity in 
relativistic plasmas does approach some terminal value, which is somewhat smaller then 
the speed of light. By determining the spectral peak of the magnetic fluctuations at each 
time step, one determined that after an initial fast growth, the typical correlation length 
of the field sets on the scalings 



(12) 



for non-relativistic and relativistic plasmas, respectively. Thus, the field scale grows 
in a self-similar manner, as one shall expect for a scale-free turbulence. We also note 
that in some respect, the field scale growth is analogous to the inverse cascade in 
two-dimensional magnetohydrodynamic (MHD) turbulence. The crucial difference is, 
however, the entirely kinetic nature of the process; at such small scales ~ c/co p the 
MHD approximation is completely inapplicable. 

RELATIVISTIC SHOCK AND ITS FORESHOCK REGION 

The simple understanding of a relativistic shock goes as follows. The particles reflected 
from the shock potential and propagate forward. They are Weibel-unstable and result 
in plasma current filamentation. The filaments interact with each other; they merge 
to form larger ones with stronger magnetic field since the current is approximately 
conserved during mergers. As filaments get stronger, they scatter particle more, thus 
causing particle heating and acceleration. At some point, the filaments are too strong, 
so they break up and result in strong longitudinal electric fields that slow down the 
ions (or the plasma in-flow, in general). Here the density builds up quickly, so this is 
the shock jump. Such a picture is, however, rather simplistic. The shock continuously 
accelerates and ejects cosmic rays (CRs) into the upstream region. Because of the non- 
linear shock acceleration, these CRs carry a substantial fraction of the shock energy, yet 
they are propagating nearly at the light speed - much faster than the shock. Hence, in 
time, the CRs deposit their energy into the magnetic field in an extremely large region 
in front. These fields will in turn affect the shock propagation and particle acceleration. 
The model of the relativistic foreshock is discussed now. 

Overall, our model of the foreshock is as follows. A shock is a source of CRs which 
move away from it, thus forming a stream of particles through the ambient medium, 
say, the interstellar medium (ISM). If the ISM magnetic fields are negligible, i.e., their 
energy density is small compared to that of CRs, the streaming instability (either the 
pure magnetostatic Weibel or the mixed-mode electromagnetic oblique Weibel-type 
instability, depending on conditions) is excited and stronger magnetic fields are quickly 
generated. The Weibel instability dispersion equation [2] has been also derived as the 
cosmic -ray streaming instability dispersion equation, see e.g., [fToh . The generated fields 
further isotropize (thermalize) the CR stream. Since less energetic particles, having a 
greater number density and carrying more energy overall, are thermalized closer to the 
shock, the generated B-field will be stronger closer to the shock and fall off away from it, 
whereas its correlation length will increase with the increasing distance from the shock. 
More energetic particles keep streaming because of their larger Larmor radii and produce 
the magnetic field further away from the shock. This process stops at distances where 
either the CR flux starts to decrease (because of the finite distance the CR particles can 
get away from a relativistic shock or because of the shock curvature causing CR density 
to decrease as r~ 2 if the shock is sub- or non-relativistic) or where the generated 
magnetic fields become comparable to the ISM field and the instability ceases. Thus, a 
large upstream region — the foreshock — is populated with magnetic fields. We now 
derive its self-similar structure. We work in the shock co-moving frame unless stated 
otherwise. 



Let's consider a relativistic shock moving along x-direction with the bulk Lorentz fac- 
tor r s h; the shock is plane-parallel and lies in the yz-plane, and x = denotes the shock 
position. The shock continuously accelerates cosmic rays, which then propagate away 
from it into the upstream region. We conventionally assume that the CR distribution over 
the particle Lorentz factor is described by a power-law: 

ncR = n (r/Yoy p+l (13) 

for 7 > 7o and zero otherwise. The index p is approximately equal to 2.2 for ultrarela- 
tivistic shocks and hq is the normalization H We assume that the above energy distribu- 
tion is the same everywhere in the upstream, that is, we neglect the nonlinear feedback 
of magnetic fields onto the particle distribution. Moreover, we assume that the forma- 
tion of the CR power-law distribution is co-temporaneous with the shock formation it- 
self, so we neglect the finite time acceleration of CRs, which may be quite important 
for higher-energy CRs. Accurate inclusion of this effect would require solution of the 
convection-diffusion equation with diffusion being calculated self-consistently from the 
self-generated fields, which are not steady at the beginning of the shock formation; all 
these issues are beyond the scope of the present paper. The CR momentum distribu- 
tion exhibits strong anisotropy: the parallel (x) components of CR momenta are much 
greater than their thermal spread in the perpendicular (yz) plane. Indeed, for a particle 
to move away from the shock, it should have the x-component of the velocity exceeding 
the shock velocity. Since both the shock and the particle move nearly at the speed of 
light, this puts a constraint on their relative angle of propagation to be less then 1 /r s h in 
the lab (observer) frame. Hence, the transverse spread of the CR particle's momenta is 
P± < P||/r s h ^ P\\ - This is also seen in numerical simulations [5|]. 

The CR particles propagate through the self-generated foreshock fields and scatter 
off them. Lower energy particles are deflected in the fields more strongly and, there- 
fore, izotropize faster than the higher energy ones, as having larger Larmor radii. At a 
position x > the CR distribution can roughly be divided into isotropic (themalized) 
component with y< y r (x) and streaming component with y> y r (x), where y r (x) is 
the minimum Lorentz factor of the streaming particles at a location x; it is also the 
maximum Lorentz factor of the randomized component at this location. The stream- 
ing component is Weibel-unstable with a very short e-folding time T ~ G^el' wnere 

1 /2 

Wp.rei = (^ne 2 n(y) / m p y) is the relativistic plasma frequency, n(y) is the density of 
streaming particles of the Lorentz factor y (tilde denotes streaming particles). Note that 
the Weibel instability growth rate depends on n of the lower density component - cosmic 
rays, in our case - measured in the center of mass frame of the streaming plasmas. For 
the lower-energy part of the CR distribution, the center of mass frame is approximately 
the shock co-moving frame, hence we evaluate the instability on the shock frame. This 
approximation is less accurate for the high-energy CR tail; however, the growth rate 

and the scale length are weak functions of the the shock Lorentz factor (<* T^ 1 ^ 2 ), so 
the result will be accurate within an order of magnitude for all reasonable values of T s h 



Conventionally the distribution is given as dn /dy°< y p with p being ~ 2.2 — 2.3 for relativistic shocks; 
hence the density of particles of energy ~ 7 is n(y) oc y~P8y<x y~P +l . 



for GRB afterglows. Here we use the proton plasma frequency because the CR electron 
Lorentz factors are about m p /m e times larger, so they behave almost like protons 
The instability is very fast: it rapidly saturates (the fields cease to grow) in a few tens of 
e-folding times x, that is in few tens of inertial lengths (also referred to as the ion skin 
length) c/co p j e i in front of the shock. Thereafter the particles keep streaming in current 
filaments and the field around them amounts to ^b ~ 0.01 — 0.001 or so of the kinetic 
energy of this group of particles: 

B 2 (r)/S7t^^ B m p c 2 rn(Y), (14) 

where is the efficiency factor obtained from PIC simulations; £b has the same mean- 
ing as the conventional Eg parameter reserved here for the ratio of the total magnetic 
energy to the total kinetic energy of the shock and which, as is seen in PIC simulations, 
is larger than near the shock because of the nonlinear evolution and filament merg- 
ers. Whereas there is a concern that the Weibel instability can be suppressed by large 
upstream magnetic fields, numerical simulations generally indicate that the upstream 
field has little influence on the Weibel instability if the field energy density does not 
exceed few percent of the total shock kinetic energy density. For our conditions with 



~ 0.01 —0.001 and reasonable CR acceleration efficiency, <^cr ~ 0.5 (see, e.g., UM) 
the instability is hardly suppressed. The correlation length of the field is of the order of 
the ion inertial (skin) length 

A(y) ~ c/co pM = (m p c 2 y/4ne 2 n(y)) 1/2 . (15) 

These random fields deflect CR particles and ultimately lead to their isotropization. 
The group of particles of the Lorentz factor y r thermalizes at the distance from the shock: 

x r ~A(y)/(2£ 5 ). (16) 

At this point, x = x r , one has y = y r by definition; no field of the strength B(y r ) and the 
scale A(y r ) can be produced at x > x r . One can show that randomization of the higher 
energy particles is small in these fields at x ~ x r , which means that these particles keep 
streaming through much larger distances x 3> x r and will produce the magnetic field 
further away from the shock. This field will be weaker and larger scale because of the 
lower density of the streaming particles h(y) <C n(y r ), according to Eqs. (fT3l)-([T5l). 

Finally, the number density of streaming CR particles at y r is n(y r ) — no(y r /yb)~ p+l . 
Therefore, 

A (y r ) ~ (m p c 2 y /A7te 2 n ) 1/2 (y,/7o) p/2 = A (y,/y ) p/2 , (17) 

where Xq is the inertial length of the lowest energy CR "plasma". Inverting this expres- 
sion yields: 

Yr ~ y [Hy r )/M 2/p ~ y (2 WAo) 2/p - (18) 

Hereafter, the subscript "r" can be omitted without loss of clarity. 

In a steady state, this field is continuously advected toward the shock (in the shock 
co-moving frame since the center of mass frame of the foreshock plasma differs from the 




FIGURE 3. A schematic representation of the foreshock magnetic fields: the coherence length is 
increasing with the upstream distance. Below are schematic graphs showing variation of the spectrum 
of the streaming part of cosmic rays and the corresponding self-generated fields (highlighted). 



shock frame) and may affect the onset and the saturation level of the Weibel instability. 
In addition, the current filaments producing the fields merge with time, so that B and A 
change while being advected. These nonlinear feed-back effects are difficult to properly 
account for in a theoretical model; hence they are omitted in the current study. PIC 
simulations can help us to quantify the effects as well as to confirm or disprove our 
assumption that the shock and the foreshock do form a self-sustained, steady state 
structure. 

The self-similar structure of the foreshock [12] immediately follows from Eqs. (PT3T) . 
CO}, (fT6l) and (fT8l) . The magnetic field correlation length is proportional to the upstream 
distance from the shock, 

X{x)~x{2$ B ), (19) 
and its strength decreases with the distance as 

B(x)~B (x/x r ip - 2)/p , (20) 

1 /2 1 /2 

where B = (Sn^ B m p c 2 n y ) and x = Xq/{2E, b ) = (m p c 2 yo/4ne 2 n ) /(2£b). In 
this estimate we neglected the advected fields B(y) as sub-dominant compared to B(y r ) 
for 7 > y r . We note here that the idea of self- similarity of the Weibel turbulence has been 



first proposed by Medvedev et al. [9] and then further elaborated by Katz et al. lfl3ll . 
whose results are in agreement with the above scalings. The eg parameter expresses 
the field energy normalized to the shock kinetic energy. The energy of cosmic rays is 
Uqr = J n(y/ 7o) (m p c 2 Y) d(y/ 7o) ~ m p c 2 Yono and constitutes a fraction <^cr of the total 
shock energy, U^. The efficiency of cosmic ray acceleration, £cr, can be as high as 
several tens percent, perhaps, up to <^cr ~ 0.5, as follows from the nonlinear shock 
modeling llllll . The scaling of £g is: 

ft~&R& (x/x )- 2{p - 2)/p . (21) 

These scalings hold while the shock can be treated as planar and while the ISM 
magnetic fields are negligible compared to the Weibel-generated fields. If the shock 
is relativistic, CR particles can occupy a narrow region in front of it. Assuming CR to 
propagate nearly at the speed of light, their front is ahead of she shock at the distance 
Ar' ~ Ar/r s h ~ R/ (2r s h) in the shock frame. Also, when the radial distance in the lab 
frame Ar = x/T s h becomes comparable to the shock radius Ar ~ R the curvature of the 
shock can no longer be neglected: the density of CR particles, which was assumed to 
be constant in our model, starts to fall as oc r 2 . This leads to a steeper decline of B 
with distance. Obviously, the first constraint is more stringent for a relativistic shock, 
whereas both are very similar (within a factor of two) for a non-relativistic shock. Hence 
one should use the first constraint. Here we also assumed that the ambient magnetic field 
is zero. It is a good approximation for usual ISM conditions II 1211 . 

We can also estimate the magnetic field spectrum at and after the shock jump as 
long as dissipation is not playing a role. The magnetic field of different correlation 
scales created in the foreshock is advected toward the shock, so a broad spectrum is 
accumulated: 

5 A ocA-^- 2 ^~A- - 091 , (22) 

where Eqs. (fT9l) and (1201) were used and p ~ 2.2 was assumed. 

We want to note that a number of simplifying assumptions has been made in our 
analysis. In particular, nonlinear feedback effects of the upstream magnetic field on the 
particle distribution, on the shock structure and on Fermi acceleration were omitted. 
The inclusion of these effects is hardly possible in any analytical model. We neglected 
in this model that the accelerated electrons can loose their energy and cool down in the 
self-generated fields while the protons keep their energy, thus causing electric fields and, 
possibly, currents to build up in the region and modify the foreshock structure. We also 
assumed that a steady state exists for the shock-foreshock system at hand. Apparently, it 
is not at all clear whether the steady state is at all possible or the system exhibits an inter- 
mittent behavior. One can envision a scenario in which the CRs overproduce upstream 
magnetic fields leading to enhanced particle scattering and the overall preheating of the 
ambient medium, which, in turn, can cause the shock to weaken, disappear and then 
re-appear in a different place further upstream. Presently available 2D PIC simulations 
of an electron-position shock do show the upstream field amplification and no steady 
state has been achieved: both upstream and downstream fields continue to grow for the 
duration of the simulations though no shock re-formation has so far been observed. We 
argue that extensive PIC or/and hybrid simulations of a shock are imperative for further 
study. 



CONCLUSIONS 



What is more important: a shock or its foreshock? Conventionally, one believes that the 
radiation observed from gamma-ray bursts is produced in the shock downstream. Indeed, 
the magnetic fields there are present, though the small-scale ones decay relatively fast. 
The foreshock fields advected into the downstream may increase the meso- scale fields 
which are not coupled to dissipation efficiently. The shock front itself can produce 
substantial radiation for certain parameters of a GRB explosion. However, the fields at 
the shock are mostly those created in the foreshock too and then advected with the flow. 
The foreshock region itself is very large in size ~ in the observer's frame and is 

filled with the magnetic field of sub-Gauss strength (a fraction of the CR energy density) 
and of correlation scale as large as ~ R/^b ~ 0.017? ~ 10 16 cm for R ~ 10 18 cm. The 
present model is very simple and neglects a number of important nonlinear feedback 
effects. Extensive simulations are needed to accurately study the foreshock structure. 
How much radiation is emitted from the foreshock region strongly depends on these 
details. It is premature to draw any conclusions on this matter. Summarizing all the 
above, we believe that the study (e.g., using PIC and other simulations) of a foreshock 
is of crucial importance for our understanding of observational signatures of relativistic 
shocks and it shall be given the top priority. 
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